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What is claimed is: 

1 . A method for reconstructing an image of a scattering medium, comprising: 
directing energy into the scattering medium at a source location on the 
scattering medium; 



predicting the energy emerging from the scattering medium using an 
equation of radiative transfer, wherein the prediction is a function of the initial guess; 

generating an objective function based on a comparison of the prediction 
with the measurement; 

generating a gradient of the objective function by a method of adjoint 

differentiation; 

modifying the initial guess of the properties of the scattering medium 
based on the gradient of the objective function; and 

generating an image representation of the internal properties of the 
scattering medium. 

2. The method according to claim 1, further comprising repeating the 
predicting of the energy emerging from the scattering medium based on the modified 
initial guess, generating the objective function and modifying the initial guess, until at 
least one of a predetermined number of repetitions has occurred and the objective 
function reaches a predetermined threshold. 



measuring the energy emerging from the scattering medium at a detector 



location on the scattering medium; 



selecting an initial guess of internal properties of the scattering medium; 
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3. The method according to claim 1, wherein the prediction depends on the 
boundary conditions. 

4. The method according to claim 3, wherein the boundary conditions 
account for a refractive mismatch at an interface between the medium and at least one of 
the detectors and source. 

5. The method according to claim 1, wherein the prediction comprises an 
iterative process producing intermediate results. 

6. The method according to claim 5, wherein the intermediate results of the 
prediction are stored. 

7. The method according to claim 6, wherein generating the gradient of the 
objective function by adjoint differences uses the intermediate results of the prediction. 

8. The method according to claim 7, wherein generating the gradient 
comprises stepping backward through the intermediate results of the prediction. 

9. The method according to claim 1, wherein the equation of radiative 
transfer is time independent. 

10. The method according to claim 9, wherein the time independent equation 
of radiative transfer is: 
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3 coW(r J co) + {\x a + ^)*F(r,co) = S(r,o) + ^ ^^co^'^Cr^)^' 

4 where ^(r,^) is the radiance at the spatial position r directed into a unit 

5 solid angle <s>, S(r,w) is the energy directed into the medium at spatial position r into a 

6 unit solid angle a, /u s is the scattering coefficient, is the absorption coefficient and 

7 pico^af) is the scattering phase function. 

1 11. The method according to claim 10, wherein the scattering phase function 

2 is: 

n 3 p(COS6) = — T— —371 

~ 2(1 + # -2gcosG) 

PI 4 where #is the angle between the two unit solid angles a> and co ', and g is 

111 5 the anisotropy factor. 

W 



1 12. The method according to claim 1, wherein the equation of radiative 

2 transfer is time dependent. 

1 13. The method according to claim 12, wherein the time dependent equation 

2 of radiative transfer is: 

3 1 ^OvM = S (r,<»,t)-a> . V¥(r,©,/) - <ti a + M f) + \i 5 g n p(fl),fl)W^V)A)' 

4 where ^r^t) is the radiance at the spatial position r directed into a unit 

5 solid angle a>, S(r,w, t) is the energy directed into the medium at spatial position r into a 
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unit solid angle co, /u s is the scattering coefficient, // a is the absorption coefficient and 
p{o),co t ) is the scattering phase function. 

14. The method according to claim 13, wherein the scattering phase function 

is: 



the anisotropy factor. 

15. The method according to claim 1, wherein the properties include at least 
one of a scattering coefficient, an absorption coefficient, an anisotropy factor, and a 
scattering phase function. 

16. The method according to claim 1, wherein the objective function is a 
normalized comparison of the predicted energy and the measured energy 

17. The method according to claim 1, wherein the objective function is based 
on the normalized sum of the differences between the predicted energy and the measured 
energy for each source detector pair, wherein a source detector pair is formed between 
each source location and each detector location. 



/?(cos0) = 



2(l + £ 2 -2gcos6) : 



t 3/2 



where #is the angle between the two unit solid angles to and co and g is 



18. The method according to claim 1, wherein the objective function is: 
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where Mi represents the actual measurements and the Pi represents the 



predicted measurements for each source defector pair /, m is the number of source 
detector pairs, where a source detector pairs is formed between each source location and 
each detector location. 

19. The method according to claim 1, further comprising minimizing the 
objective function. 

20. The method according to claim 19, wherein minimizing the objective 
function includes a one dimensional line search. 

21 . The method according to claim 20, wherein the one dimensional line 
search is performed along a direction of the gradient of the objective function. 

22. The method according to claim 20, wherein the one dimensional line 
search is performed along a gradient-dependent direction. 

23. The method according to claim 1, wherein the energy comprises near 
infra-red energy. 
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1 24. The method according to claim 1, wherein the scattering medium contains 

2 regions wherein the scattering coefficients are not substantially greater than the 

3 absorption coefficients. 

1 25. The method according to claim 1, wherein the scattering medium contains 

2 a low scattering region embedded in a high scattering region. 

1 26. The method according to claim 1, wherein the predicted energy is 

2 determined using finite element methods. 



^ 1 27. The method according to claim 1 , wherein the predicted energy is 
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2 determined using finite difference methods. 



1 28. A method for imaging the spatial optical properties of tissue, comprising: 

ill 2 (a) directing energy into the scattering medium at a source location on 

^ 3 thetiss ue; ^ 

4 (b) measuring the energy emerging from the scattering medium at a 

5 detector location on the tissue; 

6 (c) selecting and initial guess of the spatial optical properties of the 

7 tissue; 

8 (d) predicting the energy emerging from the tissue using an equation 

9 of radiative transfer in an iterative process, wherein the prediction is a function of the 
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10 initial guess and a refraction index mismatch at a boundary of the tissue, and the iterative 

1 1 process generates a plurality of intermediate predictions; 

12 (e) generating an objective function based on a normalized 

13 comparison of the prediction with the measured energy emerging from the scattering 

14 medium; 

15 (f) generating a gradient of the objective function by adjoint 

1 6 differentiation; 

1 7 (g) modifying the initial guess of the spatial properties of the tissue 

18 based on the gradient of the objective function; 

0 19 (h) repeating steps (d) through (g) until at least one of a threshold of 
N 20 modifications to the initial guess is reached and the objective function reaches a 

21 threshold; and 

1 . | 

% j 22 (j) generating an image representation of the spatial optical properties 

23 of the tissue. 



h 

ill 
I1J 

1 29. A system for reconstructing an image of a scattering medium, comprising: 



2 a source for directing energy into the scattering medium at source location on the 

3 scattering medium; 

4 a detector for measuring the energy emerging from the scattering medium at a 

5 detector location on the scattering medium; 

6 an initial guess of internal properties of the scattering medium; 

7 means for predicting the energy emerging from the scattering medium using an 

8 equation of radiative transfer, wherein the prediction is a function of the initial guess; 
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means for generating an objective function based on a comparison of the 
prediction with the measurement; 

means for generating a gradient of the objective function by a method of adjoint 
differentiation; 

means for modifying the initial guess of the properties of the scattering medium 
based on the gradient of the objective function; and 

means for generating an image representation of the internal properties of the 
scattering medium. 

30. The system according to claim 1, further comprising means for repeating 
the predicting of the energy emerging from the scattering medium based on the modified 
initial guess, generating the objective function and modifying the initial guess, until at 
least one of a predetermined number of repetitions has occurred and the objective 
function reaches a predetermined threshold. 

3 1 . The system according to claim 1 , wherein the prediction depends on the 
boundary conditions. 

32. The system according to claim 3 1 , wherein the boundary conditions 
account for a refractive mismatch at an interface between the medium and at least one of 
the detectors and source. 
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33. The system according to claim 1, wherein the prediction comprises an 
iterative process producing intermediate results. 

34. The system according to claim 33, wherein the intermediate results of the 
prediction are stored. 

35. The system according to claim 34, wherein generating the gradient of the 
objective function by adjoint differences uses the intermediate results of the prediction. 

36. The system according to claim 35, wherein generating the gradient 
comprises stepping backward through the intermediate results of the prediction. 

37. The system according to claim 1, wherein the equation of radiative 
transfer is time independent. 

38. The system according to claim 37, wherein the time independent equation 
of radiative transfer is: 



where is the radiance at the spatial position r directed into a unit solid 

angle co, S(r,w) is the energy directed into the medium at spatial position r into a unit 
solid angle co, /j s is the scattering coefficient, fu a is the absorption coefficient and p(a? 9 a? r ) 
is the scattering phase function. 



(0 W(r,co) + (ii a + \i s y¥(r,a>) = S(r,co) + \i s jjf p(p 9 v')V(r 9 rt)da> • 
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1 39. The system according to claim 38, wherein the scattering phase function 

2 is: 

3 p(cos9) = ; 777 

FK ' 2(l + g 2 -2gcos0) 3/2 

4 where #is the angle between the two unit solid angles co and co ', and g is the 

5 anisotropy factor. 

1 40. The system according to claim 1, wherein the equation of radiative 

2 transfer is time dependent. 



1 41. The system according to claim 40, wherein the time dependent equation of 

2 radiative transfer is: 



01 

I1J 

jd! 3 1 3F(gD,0 = 5(/%t0;0 _ a . VT(r,co,/)- (M a +^,)^co,/) +MJo n Kco,o)') v I , ('' J G>',/)^(B' 

i» 4 where *F(r,o,t) is the radiance at the spatial position r directed into a unit solid 

a 

:= 1 5 angle *y, S(r,w, t) is the energy directed into the medium at spatial position r into a unit 

ill 

l=| 6 solid angle co, fi s is the scattering coefficient, /i a is the absorption coefficient and p(co,o?) 

i=fc 

7 is the scattering phase function. 



1 42. The system according to claim 41, wherein the scattering phase function 

2 is: 

3 p(cos9) = - y a —777 

FK } 2(l + g 2 -2gcos6) 3/2 
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where #is the angle between the two unit solid angles co and co and g is the 
anisotropy factor. 

43. The system according to claim 1, wherein the properties include at least 
one of a scattering coefficient, an absorption coefficient, an anisotropy factor, and a 
scattering phase function. 

44. The system according to claim 1, wherein the objective function is a 
normalized comparison of the predicted energy and the measured energy. 

45. The system according to claim 1, wherein the objective function is based 
on the normalized sum of the differences between the predicted energy and the measured 
energy for each source detector pair, wherein a source detector pair is formed between 
each source location and each detector location. 

46. The system according to claim 1, wherein the objective function is: 



where M\ represents the actual measurements and Pj represents the predicted 
measurements for each source detector pair, m is the number of source detector pairs, 
where a source detector pairs is formed between each source location and each detector 
location. 



1 m ? 
cp=-Z(^-M / ) 2 
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47. The system according to claim 1, further comprising minimizing the 
objective function. 

48. The system according to claim 47, wherein minimizing the objective 
function includes a one dimensional line search. 

49. The system according to claim 48, wherein the one dimensional line 
search is performed along a direction of the gradient of the objective function. 

50. The system according to claim 49, wherein the one dimensional line 
search is performed along a gradient-dependent direction. 

51. The system according to claim 50, wherein the energy comprises near 
infra-red energy. 

52. The system according to claim 1, wherein the scattering medium contains 
regions wherein the scattering coefficients are not substantially greater than the 
absorption coefficients. 

53. The system according to claim 1, wherein the scattering medium contains 
a low scattering region embedded in a high scattering region. 
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54. The system according to claim 1, wherein the predicted energy is 
determined using finite element methods. 

55. The system according to claim 1, wherein the predicted energy is 
determined using finite difference methods. 

56. A system for imaging the spatial distribution of optical properties of 
tissue, comprising: 

(a) a source for directing energy into the scattering medium at a source 
location on the tissue; 

(b) a detector for measuring the energy emerging from the scattering medium 
at a detector location on the tissue; 

(c) an initial guess of spatial optical properties of the tissue; 

(d) means for predicting the energy emerging from the tissue using an 
equation of radiative transfer in an iterative process, wherein the prediction is a function 
of the initial guess and a refraction index mismatch at a boundary of the tissue, and the 
iterative process generates a plurality of intermediate predictions; 

(e) means for generating an objective function based on a normalized 
comparison of the prediction with the measured energy emerging from the scattering 
medium; 

(f) means for generating a gradient of the objective function by adjoint 
differentiation; 
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17 (g) means for modifying the initial guess of the spatial properties of the tissue 

1 8 based on the gradient of the objective function; 

19 (h) means for repeating steps (d) through (g) until at least one of a threshold 

20 of modifications to the initial guess is reached and the objective function reaches a 

21 threshold; and 

22 (j) means for generating an image representation of the spatial optical 

23 properties of the tissue. 

1 57. Computer executable software code stored on a computer readable 

CI 2 medium, the code for reconstructing an image of a scattering medium, comprising: 
N 3 code to direct energy into the scattering medium at a source location on 



111 



4 the scattering medium; 

5 code to measure the energy emerging from the scattering medium at a 

6 detector location on the scattering medium; 
[-»■ 

111 7 code to receive an initial guess of internal properties of the scattering 

111 

Cj 8 medium; 

9 code to predict the energy emerging from the scattering medium using an 

10 equation of radiative transfer, wherein the prediction is a function of the initial guess; 

1 1 code to generate an objective function based on a comparison of the 

12 prediction with the measurement; 

13 code to generate a gradient of the objective function by a method of 

14 adjoint differentiation; 
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15 code to modify the initial guess of the properties of the scattering medium 

16 based on the gradient of the objective function; and 

17 code to generate an image representation of the internal properties of the 

18 scattering medium. 

1 58. Computer executable software code stored on a computer readable 

2 medium, the code for imaging the spatial distribution of optical properties of tissue, 

3 comprising: 

4 (a) code to direct energy into the scattering medium at a source location on the 



O 5 tissue; 



6 (b) code to measure the energy emerging from the scattering medium at a detector 



7 location on the tissue; 
%l 8 (c) code to receive an initial guess of spatial optical properties of the tissue; 

IE 

£j 9 (d) code to predict the energy emerging from the tissue using an equation of 

k- 

111 10 radiative transfer in an iterative process, wherein the prediction is a function of the initial 

ijj 

CJ 1 1 guess and a refraction index mismatch at a boundary of the tissue, and the iterative 

12 process generates a plurality of intermediate predictions; 

13 (e) code to generate an objective function based on a normalized comparison of 

14 the prediction with the measured energy emerging from the scattering medium; 

15 (f) code to generate a gradient of the objective function by adjoint differentiation; 

16 (g) code to modify the initial guess of the spatial properties of the tissue based on 

17 the gradient of the objective function; 
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18 (h) code to repeat steps (d) through (g) until at least one of a threshold of 

19 modifications to the initial guess is reached and the objective function reaches a 

20 threshold; and 

21 (j) code to generate an image representation of the spatial optical properties of the 

22 tissue. 

1 59. A computer readable medium having computer executable software code 

2 stored thereon, the code for reconstructing an image of a scattering medium, comprising: 

3 code to direct energy into the scattering medium at a source location on 
Q 4 the scattering medium; 

5 code to measure the energy emerging from the scattering medium at a 

Q\ 

\l 



6 detector location on the scattering medium; 
'f\ 7 code to receive an initial guess of internal properties of the scattering 

| 8 medium; 

9 code to predict the energy emerging from the scattering medium using an 

10 equation of radiative transfer, wherein the prediction is a function of the initial guess; 

1 1 code to generate an objective function based on a comparison of the 

12 prediction with the measurement; 

13 code to generate a gradient of the objective function by a method of 

14 adjoint differentiation; 

15 code to modify the initial guess of the properties of the scattering medium 

16 based on the gradient of the objective function; and 
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code to generate an image representation of the internal properties of the 



scattering medium. 

60. A computer readable medium having computer executable software code 
stored thereon, the code for imaging the spatial distribution of optical properties of tissue, 
comprising: 

(a) code to direct energy into the scattering medium at a source location on the 

tissue; 

(b) code to measure the energy emerging from the scattering medium at a detector 
location on the tissue; 

(c) code to receive an initial guess of spatial optical properties of the tissue; 

(d) code to predict the energy emerging from the tissue using an equation of 
radiative transfer in an iterative process, wherein the prediction is a function of the initial 
guess and a refraction index mismatch at a boundary of the tissue, and the iterative 
process generates a plurality of intermediate predictions; 

(e) code to generate an objective function based on a normalized comparison of 
the prediction with the measured energy emerging from the scattering medium; 

(f) code to generate a gradient of the objective function by adjoint differentiation; 

(g) code to modify the initial guess of the spatial properties of the tissue based on 
the gradient of the objective function; 

(h) code to repeat steps (d) through (g) until at least one of a threshold of 
modifications to the initial guess is reached and the objective function reaches a 
threshold; and 
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21 (j) c °de to generate an image representation of the spatial optical properties of the 

22 tissue. 
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